function [v1, p_exit] = valfunc_I(v2, s,y,param,glob,options)
% inputs 

    F            = menufun('F',s,y,param,glob,options);
              
    L            = menufun('L',s,y,param,glob,options);   
    
    adj_cost     = menufun('C',s,y,param,glob,options);
    
    sa           = [L, s(:,2)];
    
    cv           = param.sdf*(1 - param.gamma)*v2(sa); 

    
    p_exit       = 1 - logncdf(cv, param.mu_f, param.sigma_f);
    idx          = p_exit < 1;
    
    exp_f_cost   = zeros(size(p_exit));
    exp_f_cost(idx)   = cond_logn_exp(cv(idx), param.mu_f, param.sigma_f);
    
    v1           = F - adj_cost + (1-p_exit).*(cv - exp_f_cost);

    
end